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\0 ■ ABSTRACT 

o ■ 
o : 

CN . We compare a large set of cosmologies with WMAP data, performing a fit based 

O ! on a MCMC algorithm. Besides of ACDM models, we take dynamical DE models, 

Q , where DE and DM are uncoupled or coupled, both in the case of constant coupling 

I and in the case when coupling varies with suitable laws. DE however arises from a 

^Sj ' scalar field self-interacting through a SUGRA potential. We find that the best fitting 
model is SUGRA dynamical DE, almost indipendently from the exponent a in the 

^ ■ self-interaction potential. 

. The main target of this work are however coupled DE models, for which we find 

\^ , limits on the DE-DM coupling strength. In the case of variable coupling, we also find 

I that greater values of the Hubble constant are preferred. 

^ : 

■ Subject headings: cosmology: theory - dark energy 

Pi 1. Introduction 

■ Cosmologies with flo,de — 0.75, ^lo,m — 0.25, ^lo,b — 0.03, h ~ 0.7 and ~ 1 fit available data 
J> . on CMB, LSS, and SNIa. Respectively, the above symbols are: the present density parameters 

^ , for dark energy (DE), whole non-relativistic matter, baryons, the Hubble parameter in units of 

^ [ 100 km/s/Mpc, and the primeval spectral index. CMB is the Cosmic Microwave Background, LSS 

is the Large Scale Structure. Possible references for the above data are: Tegmark et al. (2001), 
De Bernardis et al. (2000), Hanany et al. (2000), Halverson et al. (2002), Spergel et al. (2003), 
Percival et al. (2002), Efstathiou et al. (2002), Riess et al. (1998), Perlmutter et al. (1999). 

The nature of DE is not clear. This paper is devoted to comparing various hypotheses on its 
nature among them and with WMAP data. It is already known that, in order to fit data, DE must 
have a largely negative pressure: pde = wpde-, with —w ~ 0.7-1 {pde is DE energy density). False 
vacua obey such a state equation with w = —1, but assuming DE to arise from false vacuum implies 
a fine-tuning of pde, at the end of the EW phase transition by 1 : 10^^. Otherwise, DE could be 
due to a self-interacting scalar field (j) (Wetterich 1988, Ratra & Peebles 1988, RP hereafter). In 
this case we have dynamical DE (dDE) or quintessence. 



The principal analysis of WMAP first-year data (jSpergel et al. 20031 ) constrained flat ACDM 
models defined by six parameters: ^obh"^, ^omh'^, h, n^, the fluctuation amplitude A and the 
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optical depth r. As possible extensions of ACDM cosmologies, several works considered models 
with a fixed state parameter w = Pde/Pde (Spergel et al. 2003, Bean et al. 2004, Tegmark et 
al. 2004, Melchiorri et al. 2004), or adopted z-dependent paramcterizations of w{z) interpolating 
between early-time and late-time values (Corasaniti et al. 2004, Jassal et al. 2005, Rapetti et al. 
2004). A general conclusion was that current data mostly allow to constrain only the present state 
parameter, w{z = 0) < —0.80. Here we shall extend the comparison to dDE models admitting a 
tracker solution. 

Possible interactions of DE were also considered in the literature. Severe limits can be set 
to its interaction with baryons; Amendola (2000, 2003), Gasperini, Piazza &; Veneziano (2002), 
Perrotta &; Baccigalupi (2002), among others, discussed models where it interacts with Dark Matter 
(DM). Amendola & Quercellini (2003) performed a preliminary study on the limits on DE-DM 
interactions, arising from WMAP data, for DE theories with constant coupling (ccDE). It is also 
possible to consider cosmologies where the DM-DE coupling constant C, that we shall define below 
in accordance with Amendola (2000), is replaced by 1/^ and is therefore variable in time (vcDE). 
A physical origin for these models was recently discussed by Mainini Sz Bonometto (2003; see also 
Colombo et al. 2003). 

In this paper we perform a systematic fit of WMAP data to these different cosmologies: dDE, 
ccDE, vcDE, and ACDM. For the first 3 of them a SUGRA self-interaction potential will be taken 
(Brax & Martin 1999, 2000; Brax, Martin Sz Riazuelo 2000); its expression is given in Sec. 2, 
together with essential information on coupled and uncoupled DE models. Further information on 
coupled models can be found, e.g., in Colombo et al. (2004). 

The fit is performed by using a MCMC (MonteCarlo Markov Chain) algorithm, in the same 
way as the WMAP team fit the ACDM model. In Sec. 3, we describe the fitting procedure. In the 
case of that model we re-obtain their results. 

The number of parameters to be fitted depends on the nature of the model. ACDM and vcDE 
models have the same number of parameters. In dDE models there is one extra parameter and 
ccDE models have a further parameter in respect to the previous cathegory. In evaluating the 
success of a fit, however, the number of parameters is taken into account. In Sec. 4 we perform sich 
comparison and in Sec. 5 we draw some conclusions. 



2. The SUGRA potential 

All models with dynamical DE we consider make use of a SUGRA potential 

to yield the self-interaction of the scalar field cj); here rrip = is the Planck mass; V depends on 

the two parameters A and a. When they are assigned, Cl^e is however fixed. We therefore prefere 
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to use Qde and A as independent parameters. The latter scale is related to a as shown in Figure [H 
which holds also in the presence of DM-DE coupling. The potential appears in the Lagrangian 

^ = V^[\9f.udf,^d,4>-V{<p)'^ ; (2) 

in the absence of coupling, the equation of motion 

4> + 2-(j) + a^V'U) = (3) 
a 

follows, while the expressions for energy density = Pij>;kin + P^;pot and pressure = p^-kin- P<f>;pot 
are obtainable by combining the terms 

P^,kin = , P<t>,pot = V{(l3) . (4) 

In the radiation dominated era, eq. ^ admits the tracker solution 

^'^^^ = aaA'^+W , (5) 

with ga = a{a + 2)2/4(a + 6). This tracker solution fairly defines the initial conditions also in the 
presence of coupling, although, in the presence of coupling, the equation of motion becomes 

(7^ + 2-0 + aV(0) = CaVc (6) 
a 

while CDM energy density varies according to the equation 

Pc + 3-pv = -Ca^pc ; (7) 
a 

baryon and radiation equations keep the usual shape. The coefficient C sets the coupling strength 
and, in the case of constant coupling (ccDE), we can set 

so to parametrize the coupling strength through the adimensional coefficient /3. Here we consider 
also a (/)-dependent coupling (vcDE) 

c = (9) 

for which there is no extra parameter to fix the coupling strength. 

3. Comparison with WMAP data through a MCMC algorithm 

WMAP data have been extensively used to constrain cosmological parameters. They are high 
precision estimates of the anisotropy power spectrum Cf^ up to / ~ 900, as well of the TE correlation 
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power spectrum Cf^ up to Z ~ 450. We shall fit these data to possible cosmologies, by considering 
a parameter space of 6 to 8 dimensions. A grid-based likelihood analysis would require a huge CPU 
time and we use a Markov Chain Monte Carlo (MCMC) approach, as it has become customary for 
CMB analysis (Christensen et al. 2001, Knox et al. 2001, Lewis et al. 2002, Kosowski et al. 2002, 
Dunkley et al. 2004). 

To this aim, as in any attempt to fit CMB data to model s, a linear code providing C/'s is 



needed. Here we use our optimized extension of CMBFAST (jSeljak &: Zaldarriaga 19961 ). able 



to inspect also ccDE and vcDE cosmologies. The likelih ood of each mod el is then evaluated 



throu gh the publicly available code by the W MAP team (j Verde et al. 20031 ) and accompanying 



data (jHinshaw et al. 20031 : iKogut et al. 20031 ). 



A MCMC algorithm samples a known distribution /2(x) by means of an arbitrary trial distri- 
bution p(x). Here £ is a likelihood and x is a point in the parameter space. The chain is started 
from a random position x and moves to a new position x', according to the trial distribution. The 
probability of accepting the new point is given by £(x')/£(x); if the new point is accepted, it is 
added to the chain and used as the starting position for a new step. If x' is rejected, a replica of x 
is added to the chain and a new x' is tested. 

In the limit of infinitely long chains, the distribution of points sampled by a MCMC describes 
the underlying statistical process. Real chains, however, are finite and convergence criteria are 
critical. Moreover, a chain must be required to fully explore the high probability region in the 
parameter space. Statistical properties estimated using a chain which has yet to achieve good con- 
vergence or mixing may be misleading. Several methods exist to diagnose mixing and convergence, 
involving either single long chains or multiple chains starting from well separated points in the 
parameter space, as the one used here. Once a chain passes convergence tests, it is an accurate 
representation of the underlying distribution. 

In order to ensure mixing, we run six chains of ~ 30000 points each, for each model cate- 
gory. We diagnose convergence by requiring that, for each parameter, the variances both of the 
singl e chains and of the whole set of chains {W and B, respectively) satisfy the Gelman & Rubin 



test (jVerde et al. 20031 : ICelman Rubin 19921 ). R < 1.1 with: 



[{N-l)/N]W + il + l/M)B 

R = ^ ■ (10) 

Here each chain has 2A^ points, but only the last N points are used to estimate variances, and 
M is the total number of chains. In most model categories considered, we find that the slowest 
parameter to converge is A. 



4. Results 



The basic results of this work are shown in the Tables [TH3l For each model category we provide 
the expectation values and the variance of each parameter, as well as the values of the parameters 
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of the best fitting models. The corresponding marginahzed distributions are plotted in figures [2H11 
while joint 2D confidence regions are shown in figures [5HZ1 

The values of can also be compared, taking into account the number of degrees of freedom. 
This is shown in Table [H The smallest is obtained for the uncoupled SUGRA model, which 
performs slightly better than ACDM. 

Tables 1 and 2 and the corresponding figures, concerning uncoupled or constant-coupling 
SUGRA models, show that WMAP data scarcely constrain A = log(A/GeV), allowed to vary 
from ~ — 12 to 16. On the contrary, in the case of a coupling, loose but precise limitations on 
A are set, as is shown in Table 3. In the presence of this coupling, the 2-a A-interval ranges from 
~ 10 to ~ 3 • 10^°GeV. 

These constraints apparently arise from polarization data: when DE begins to affect cosmic 
expansion, the cp evolution acts on DE density and pressure, as well as on its coupling, and requires 
an /-dependence of multipole amplitudes ablt to meet data only for a limited range of A values. 
This can be noticed in Figures [8] and [H where the best-fit and Cf^^ spectra for all best-fit 
models (apart of ACDM) are compared. 

Fig. [8] also shows why no model neatly prevails. At large I all best-fit models yield similar 
behaviors. Discrimination could be improved only through better large angular scale observations, 
especially for polarization, able to reduce errors on small-/ harmonics. 



5. Discussion 

5.1. Uncoupled SUGRA models 

SUGRA uncoupled models are found to be consistent with WMAP data. The ratio w = p/p, 
for these models is typically < —0.80 at z = 0. However, they exhibit a fast variation of w, which 
is already ~ —0.6 at z ~ 1-2. Such a decrease, however, does not cause a conflict with data and 
these models perform even better than ACDM. 

We also considered the effect of adding some priors. For uncoupled or constant-coupling 
SUGRA models, the analysis in the presence of priors leads to analogous conclusions. There are 
however variations in the estimate of cosmological parameters. First, the opacity r is pushed to 
values exceeding the ACDM estimates (see also Corasaniti at al. 2004). This can be understood 
in two complementary ways: (i) ISW effect is stronger for dDE than for ACDM, as the field (p 
evolves during the expans ion; hence, DE effects extend to a greater redshift, increasing Cf values 



in the low-/ plateau (e.g., IWeller &: Lewis 20031 ). To compensate this effect the fit tends to shift Us 
greater values. Owing to the r-n degeneracy, this is then balanced by increasing r. (ii) It must also 
be reminded that the expected TE correlation for dDE, at low /, is smaller (Colombo et al. 2003). 
Then, when TE correlation data is included, a given correlation level is fit by a greater r. In any 
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Table 1: SUGRA parameters for dDE models: for each parameter the expectation value (x), its 
variance CTx, and its maximum likelihood value Xmax, in the 7-dimensional parameter space, are 
shown. 



X 


{x) 




•^max 




0.025 


0.001 


0.026 


O ?)2 


0.12 


0.02 


0.11 


h 


0.63 


0.06 


0.58 


T 


0.21 


0.07 


0.28 


ris 


1.04 


0.04 


1.08 


A 


0.97 


0.13 


1.11 


A 


3.0 


7.7 


13.7 



Table 2: SUGRA parameters for ccDE models; the parameter space is 7-dimensional and parameter 
are shown as in the previous Table. 



X 


{x) 


C^x 


^max 




0.024 


0.001 


0.024 


O ^2 


0.11 


0.02 


0.12 


h 


0.74 


0.11 


0.57 


T 


0.18 


0.07 


0.17 


Us 


1.03 


0.04 


1.02 


A 


0.92 


0.14 


0.93 


A 


-0.5 


7.6 


8.3 


3 


0.10 


0.07 


0.07 



Table 3: SUGRA parameters for vcDE models. At variance from other model categories, A here is 
constrained. Parameter values are shown as in Table 1. 



X 


{x) 


(^x 


^max 




0.025 


0.001 


0.026 


^''o,dm''' 


0.11 


0.02 


0.09 


h 


0.93 


0.05 


0.98 


T 


0.26 


0.04 


0.29 


n. 


1.23 


0.04 


1.23 


A 


1.17 


0.10 


1.20 


A 


4.8 


2.4 


5.7 
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2 4 6 8 10 12 

Iog(A/GeV) 

Fig. 1. — Values of a corresponding to given A scales in SUGRA models with ^o^dm = 0.27. The 
tiny differences, just above noise level, for coupled models {ap and aij^ for constant coupling and 
(f)~^ coupling, respectively), are also shown. 



Table 4: Goodness of fit. Wc lists the number of degrees of freedom (d.o.f.), the reduced Xe//i ^'^d 
the corresponding probability of the best-fit model. ACDM models have 1342 degrees of freedom, 
uncoupled and vcDE models have 1341, while ccDE models have 1340. 







prob. 


dDE 


1.062 


5.2 % 


ccDE 


1.066 


4.7% 


vcDE 


1.074 


2.9 % 


ACDM 


1.066 


4.7% 
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Fig. 2. — Marginalized distributions for the 7-parameters SUGRA model with no priors (solid 
lines), BBNS prior (long dashed) or HST prior (dot-dashed). Short dashed (dotted) vertical lines 
show the boundaries of 68.3 % c.l. (95.4 % c.l.) interval; for A only upper limits are shown. 
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Fig. 3. — As Fig. [2] but for the 8-parameters ccDE model. For A and (3 only the upper c.l. 
boundaries are shown. 
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Fig. 4. — As Fig. [2] but for the 7-parameters vcDE model. 
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Fig. 5. — Joint 2D constraints for SUGRA models. Light (dark) shaded areas delimit the region 
enclosing 68.3 % (95.4 %) of the total points. 
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Fig. 6. — As Fig. [5] but for ccDE models. 
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Fig. 7. — As Fig. [5] but for vcDE models. Here, cosmological parameters are more stringently 
constrained than in other models. 
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1 10 100 1000 

1 

Fig. 8. — Cf spectra for the best fit SUGRA (solid line), ccDE (dotted line) and vcDE (dashed). 
Dual-axion (dot-dashed) models are also shown; they can be considered a particular case within 
vcDE models, for which 10 < A < 10.5 (see text). The binned first-year WMAP data are also 
plotted. 
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Fig. 9.— Best fit Cf ^ spectra. 
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Fig. 10.— Best fit Cf spectra. 
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case, however, r ~ 0.07 keeps consistent with data within less than 2-(j's. 

Greater r's affect also ri^/i^ estimates, whose best-fit value is then greater, although con- 
sistent with ACDM w ithin 1-a. Adding a prior on Q^h^ = 0.0214 ± 0.0020 (BBNS estimates 



(jKirkman et al. 20031 )) lowers h, slightly below HST findings, but still well within 1-a. We there- 



fore consider the effect of adding a prior also on h. 

The effects of these priors are shown in Figures [2] and [3l The new distributions are given by the 
dashed red line (prior on Qft/i^) and the dot-dashed blue line (prior on h). The former prior affects 
mainly r and n^; t and Ug are lowered to match WMAP's findings, the high tail of r distribution 
is partly suppressed. The physical analysis of primeval objects causing reionization (Ciardi et al. 
2003, Ricotti et al. 2004). which are still allowed but certainly not required. 



5.2. ccDE SUGRA models 

The latter prior favors greater h values. In the absence of coupling, this favors low-A models, 
closer to ACDM. In fact, the sound horizon at decoupling is not affected by the energy scale A, 
while the comoving distance to last scattering band is smaller for greater A's. Then, for greater 
A's lower h values are favored, so to meet the first peak. In the presence of coupling, there is a 
simultaneous effect on /?: greater f5's yield a smaller sound horizon at recombination, so that the 
distribution on h is smoother. 

A previous analysis of WMAP limits on constant coupling models had been carried on by 
Amendola & Quercellini (2003). Their analysis concerned potentials V fulfilling the relation 
dV/d(f) = BV^ , with suitable B and N. Furthermore, they assume that r = 0.17. Our analy- 
sis deals with a different potential and allows more general parameter variations. The constraints 
on P we find are less severe. It must be however outlined that (3 ^ 0.1-0.2 seem however forbidden 
by a non-linear analysis of structure formation (Maccio et al. 2004). 



5.3. vcDE SUGRA models 

The main peculiarity of (/)~^-models is that, although the value is not much worsened in 
respect to other cases, parameters are more strongly constrained in this caseq. This is evident in 
Fig. [71 In particular, at variance from the former case, the energy scale A is significantly constrained. 

Several other parameters are constrained, similarly to coupled or uncoupled SUGRA models. 
What is peculiar of ^-models is the range of h values which turn out to be favored. The best-fit 
2-a interval does not extend much below 0.85 . 

This is a problem for these models in their present form. Should new data support a greater 
h value, vcDE models would however be favored. Quite in general, however, we must remind that 
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vcDE models were considered here only in association with a SUGRA potential; vcDE models with 
different potentials can possibly agree with a smaller h. 

5.4. Dual axion models 

A particular but significant case, among vcDE models, are dual axion models (DAM; for a 
detailed analysis see Mainini & Bonometto 2004, Colombo et al 2004). 

These models arise from a generalization of the PQ solution of the strong— CP problem (Peccei 
& Quinn 1977), obtained by replacing the NG potential acting on a complex scalar field ^ with 
a quintessence potential admitting a tracker solution. While, in the usual PQ approach and its 
generalization, only the phase 9 of the $ field is physically sgnificant, in the DAM also its modulus 
(f) is physically observable. In both cases, quanta of the 9 field are axions and, if DM is made of 
axions, its density parameter ^o,m is set either by the FpQ parameter, in the NG potential, or by 
the energy scale A, in the quintessence potential. In the latter case, therefore, both DM and DE 
features follow from fixing the single parameter A. 

A detailed analysis of DAM shows that it belongs to the set of vcDE cosmologies; its main 
peculiarity is that A and r2„.„j arc no longer independent degrees of freedom. Henceforth, this model 
is strongly constrained and depends on the same number of parameters as SCDM. By varying ^o,m 
in its physical range, A keeps however close to 10. This great value of A therefore causes a conflict 
with the observed h range. All previous comments on vcDE models however hold, while it is also 
possible that extra contribution to DM, arising from topological singularities expected in axion 
theories, can ease the h problem. 

6. Conclusions 

The first evidences of DM are fairly old; they date some 70 years ago. However, only in the 
late Seventies limits on CMB anisotropics showed that non-baryonic DM had to be dominant. DE 
can also be dated back to Einstein's cosmological constant, although only SN la data revived it, 
soon followed by data on CMB and deep galaxy samples. 

The main topic of this paper is the fit of WMAP data with various DE models: ACDM, 
SUGRA dynamical and constant coupling DE models, as well as variable coupling DE models. In 
the last model DM and DE are coupled in a non-parametric way, with C = 4>~^. 

The fits of WMAP data to these models yield similar x^'s. At variance from other model 
categories, however, in vcDE models, CMB data constrain A. This is due to the stronger effects of 
(f) variations on the detailed ISW effect, as they affect both DE pressure and energy density, as well 
as DE-DM coupling. In principle, this strong impact of (p variation could badly disrupt the fit and 
make vcDE models significantly farther from data. It is noticeable that this does not occur. 
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The success of vcDE models would be complete if the favored range of values of the Hubble 

parameter (h ~ 0.8-1) could be slightly lowered. This range is however obtained just for a SUGRA 
potential. Furthermore, primeval fluctuations were assumed to be strictly adiabatic while, in axion 
models, a contribution from isocurvature modes can be expected. This could legitimately affect 
the apparent position of the first peak in the anisotropy spectrum, so completing the success of the 
model, in a fully self-consistent way. 
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